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S: 

Abstract. This paper is devoted to establish an invariance principle where the limit process is a multifractional 
. Gaussian process with a multifractional function which takes its values in (1/2, 1). Some properties, such as regularity 

and local self-similarity of this process are studied. Moreover the limit process is compared to the multifractional 
i—{ t Brownian motion. 

m ' 

Resume. Ce papier a pour but d'etablir un principe d'invariance dont le processus limite est gaussien et multifrac- 
tionnaire avec une fonction de Hurst a valeurs dans (1/2, 1). Des proprietes telles que la regularite et l'autosimilarite 
locale de ce processus sont etudiees. De plus, le processus limite est compare au mouvement brownien multifraction- 
naire. 
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^ ■ 1. Introduction 

X 

Fractional and multifractional processes or fields have been extensively studied because they provide relevant 
models in many situations such as mathematical finance, network traffic, physics, and other fields related to 
applied mathematics. See, e.g., [10] for a convenient reference. 

The most famous and simplest fractional process is fractional Brownian motion which has been introduced 
by Mandelbrot and Van Ness [12]. For a fixed H £ (0, 1), it can be defined as a centered Gaussian process 
Wh with covariance 



E[W H (t)W H (s)} = E[ ^ (1)2] ([tr + \s\ 2H \t s\ 2H ). 

Its main properties, which depend on H , are the following. 
• Wh is self-similar with index H . 
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• Wh admits H as Holder exponent at each point. 

• If H > 1/2 the increments of Wh (which are stationary) satisfy the long-range dependence relation as 
n — > oo 

E[W H (l)(W H (n + 1) - W H {n))] ~ H(2H - l)E[W H (l) 2 }n 2H - 2 . (1) 



Let us focus on this last property. Notice that the sequence of the increments of fractional Brownian motion 
called the fractional white noise is not the only Gaussian sequence with long-range dependence (see, for 
instance, [15], p. 336). However, the fractional white noise serves as a universal Gaussian model for long- 
range stationary phenomena. This is due to the invariance principle, which was established in [9, 15]. 

Theorem 1 (Invariance principle). Let H belong to (1/2,1), and {X n } n€ ^ be a stationary sequence of 
centered Gaussian random variables with covariance satisfying when n — > oo 

E[X X n }~cn 2H - 2 , 

with c > 0. Then defining for every t>0, 

n=l 

the finite- dimensional distributions of S N converge to those of CqWh as N goes to oo, where Wh is a 
fractional Brownian motion of index H and Cq = H^ 1 (2H — l) _1 c. 

In the previous statement, \_x\ denotes the integer part of the real number x. 

One of the main drawbacks of fractional Brownian motion for applications lies in the homogeneity of its 
properties, due to the fact that its pointwise Holder exponent is constant. Hence, multifractional processes 
and fields, for instance, the multifractional Brownian motion, have been introduced and attracted attention 
[3, 4, 6, 8, 11, 13]. Multifractional processes have locally the same properties as the fractional processes, 
but not globally. In fact their properties are not governed by a constant exponent H <= (0,1), but by a 
(0, l)-valued function h which is called the multifractional function. For instance, multifractional processes 
are locally self-similar and their pointwise Holder exponents vary along the trajectory. 

Because they take into account the variations of properties such as regularity, multifractional processes 
with a (1/2, l)-valucd multifractional function could be relevant alternatives to fractional processes with 
He (1/2,1) to provide models for long-range phenomena [2]. The main aim of this paper is to prove the 
existence of multifractional Gaussian processes which can serve as universal Gaussian models for long-range 
dependence. More precisely, we establish an invariance principle where the limit process is a multifractional 
Gaussian process with long-range increments. A study of these processes is developed. Besides local self- 
similarity and Holder regularity, a representation result is given. 

Maybe the most famous multifractional process is the multifractional Brownian motion. In this paper, we 
also compare as much as possible the limit processes we obtain with the multifractional Brownian motion. 

In Section 2, we recall the definition and some properties of the multifractional Brownian motion. The 
main result of this paper is established in Section 3. Section 4 deals with local self-similarity and regularity 
properties of the limit processes which are obtained in the main result. In Section 5, we give a representation 
of the limit processes. This representation allows us to justify a universal property for the limit processes. 
Finally, we give in Section 6 some applications of our results, in particular, an extension of results of [14]. 
The proof of some technical lemmas are postponed to the Appendices. 

Notation 

• For two random processes (or fields) Z\ and Z^, we denote by Z\ = D Z^ 1 the fact that the finite dimensional 
margins of Z\ are equal in distribution to those of Z2. 
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• For a random process (or field) Z and a sequence {Zn}n^ (rcsp. {Z e } e> o) of random processes (or fields) 
we denote by lim^^^ Zjq = Z (resp. lim^ Z £ = Z) the fact that the finite-dimensional distributions o 
Zm (rcsp. Z e ) converge to those of Z as N goes to oo (resp. e goes to 0). 



2. Some preliminaries on multifractional Gaussian processes 

Multifractional Brownian motion was the first multifractional Gaussian process introduced, independently 
in [6] and [13]. 

Let -B(d£) be a Brownian measure and B(dx) be its Fourier transform (See [15], Chapter 7 for more 
details). Multifractional Brownian motion with a multifractional function h: [0, oo) — > (0, 1) can be defined 
for every t by 

following [6]. The constant C(h(t)) is such that E[W^(1) 2 ] = 1. In this case, the covariance of the multifrac- 
tional Brownian motion is given for every t and s by [2] 

E[W h (t)W h (s)] = D(h(t), fc(«))(|t| fc W + fcW + \ s \Ht)+Hs) _ \ t _ s \h{t)+h(s)^ (3) 
where for every H\ , H 2 



_ s/Y^m + l)r(2g 2 + 1)8^(71^)8111(71^ 

1 U 2) 2T(H 1 +H 2 + l)sm(n(H 1 +H 2 )/2) 

Now we assume that h is /3-H61der continuous and 

sup h < p. (4) 

Then multifractional Brownian motion has locally the same properties as fractional Brownian motion. Mul- 
tifractional Brownian motion is locally self-similar, that is, for every t > 0, 



\ W h {t + eu)-W h {t) \ 
H { 7W) j u>0 ^ { M*) («)Uo 



where the convergence is in distribution in the space of continuous functions endowed with the topology of 
the uniform convergence on compact sets, and for every t its Holder pointwise exponent otw h (t) is almost 
surely equal to h(t). 

3. The main result 

In this paper, we consider a centered Gaussian field X = {X n (H), H G (1/2, 1), n £ N}, and two real numbers 
o and b such that 1/2 < a < b < 1. We also consider a continuous function /i : R — > [a, b] c (1/2,1). For every 
n,iVeN and t>0, we define h% = h(n/N) and 

im X (h N ) 

oN/ + \ _ \ ^ ^nV^n ) 

a h w - N h« ' 

n— 1 

Our aim is to study the asymptotic behavior of {S^}n£N as N goes to 00. We will use the following 
assumptions. 
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• Assumption (i) For every M > 0, the map 

(j,k,H 1 ,H 2 )^E[X J (H 1 )X k (H 2 )} 

is bounded on {(j, k) G N 2 , \j — k\< M} x [a, b} 2 . 

• Assumption (ii) There exists a continuous function R: [a,b] 2 — > (0,oo) such that 

lim sup \(j-k) a - ai - H *E[X j (Hi)Xk(H 2 )]-R(H 1 ,H 2 )\=0. (5) 

3-k-xyo (Hi,ff 2 )G[o,6] 2 

Assumption (ii) expresses that the field is asymptotically stationary and fractional for each H . Hence, 
for a fixed H, the random sequence nt—>X n (H) satisfies the classical invariancc principle. 

For a centered Gaussian field, X = {X n (H), H G (1/2, l),n G N} satisfying Assumption (ii); the function R 
in Assumption (ii) will be called the asymptotic covariance of X. More generally, a measurable and bounded 
function R: [a, 6] 2 — > (0, oo) such that there exists a centered Gaussian field X = {X n (H),H G (1/2, l),n G N} 
satisfying Assumption (ii) with R as asymptotic covariance will be called an asymptotic covariance. 

Fix an asymptotic covariance R. We denote by Qr the space of all centered Gaussian fields {X n (H)} n: H 
satisfying Assumptions (i) and (ii) with R as asymptotic covariance. We define for every t, s, Hi, H 2 

K(t,s;Hi,H 2 ) = R{Hi,H 2 )l t > s + R(H 2 ,Hi)l t<s , (6) 

and the function 

TV : [0, oo) 2 — > K, 

{9,a)^K(d,a;h{9),h{<j))\d-a\ h ^ +h ^- 2 . (7) 

Let us remark the symmetric property fulfilled by 1Z: for every t, s, Hi, H 2 , 

K{t,s;H l ,H 2 )=1l{s,t;H2,Hi). 

Let us also note that 72.* is locally integrable. Indeed, by the boundedness of R, on every compact set 
K C [0, oo) 2 there exists a constant c = ck, such that for every (8, a) £ K satisfying 9 ^ a 

\ll*{6,a)\<cx{\0-cj\ 2a - 2 + \e-cj\ 2h - 2 ). (8) 

Theorem 2. Under Assumptions (i) and (ii), when N goes to oo, the finite- dimensional distributions of 
Sfr converge to those of a centered Gaussian process Sh with covariance given for t, s > by: 



E[S h (t)S h (s)}= d9 dan(e,o-,h(9),h(a))\6-a\ h ^+ h ^- 2 
Jo Jo 

= [ dO [ daTZ*(6,a), (9) 
Jo Jo 

where the integral in the right-hand side of (9) is always defined. 

Proof. For every n and N, we define X n , N := X n (h%)/N h ™ and we let for every t, s > 

f dO f daK*(6,a). (10) 
o 



I(s,t) 

Because {X n (H)} n ,H is a centered Gaussian field, it is enough to show that for every t, s > 0: 

[Nt\ [Ns\ 

lim E[S£(t)S?(8)]= lim V J]E[X ktN X jtN ]=I(t,a). 

3=1 k=l 
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Formula (5) deals with the asymptotic behavior of E[Xj(Hi)Xk(H 2 )] when j — k goes to +00. Moreover, it 
also gives 

lim sup \(k-j) 2 - Sl - B 'E[X j (H 1 )X k (H 2 )]-R(H 2 ,H 1 )\=0, (11) 

so that we can write for every j , k (with j =/= k), Hi and H 2 

E[X j (H 1 )X k (H 2 )] = \j -k\ H i +H *- 2 (R,(j,k;H u H 2 ) +r {j,k;H 1 ,H 2 )), (12) 
where 

lim sup fo(J, k; Hi, H 2 ) = 0. 

|j-fc|->oo (H 1 ,H 2 )e[ Q ,6]2 

Let 77 > 0. Following (12) and the fact that a <h <b, there exists an integer M = M n such that for \ j — k\ > M 
and for every N, 

1 



(I *\ 


77 


j-k 


\N' N J 


N 2 


N 



(j_ *\ 




j-k 


\N' N J 


N 2 


N 



<E[x ktN Xj, N ] < —n 

Thanks to Lemma 2 (see the Appendix), 

[Nt\ [Ns\ 

I(s,t) - r)J{s,t) < liminf V V E[X ft>iv X,- iv]lu_ fc | >M 

3=1 fc=i 

< limsup ^ Y E [^"fc,iv^7,iv]l|j-fe|>M < J( s >*) +r)J(s,t), 



(13) 



AT— >00 



3=1 k=l 



where 



J(s,t) = / d0 f da\9-a 
Jo Jo 



h(S)+h(a)-2 



(14) 



(15) 



Now we deal with X^=i^ El=i"' ^[Xfc^-^^jvllij— fc|<Af- From Assumption (i) and a < h < b, we get for every 
j, fc and N 

\E[X k , N Xj, N }\< 1 sup lEtX^ffOX^ffa)]! 

ATi (H 1 ,JJ 2 )e[a,6]2 

^ su P |IE[^(i/i)A fe (i/ 2 )]|, 

where sup {Hi H2)e[a b]2 \E[X j (Hi)X k (H 2 )}\ <oo. Moreover, 



(16) 



LJVtJ [Ns\ 

Y l\j-k\<M<W(t + s))(l + 2M). 
3=1 fe=i 

Hence, combining (16) with (17) and using 2a > 1, we get 

[JVtJ [Ns\ 

J im E E E[^fc,JV^-,Jv]l| J -- fc |<A/ = 0. 

N— »oo z — * z — rf 

j=l fe=l 



(17) 



(18) 
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Combining with (14) and letting 77 — >• 0, wc obtain for every t and s, 
lim K[S?(t)S£(s)]=I(t,s). 

iV— >oo 

Hence, the finite-dimensional distributions of converge to those of a centered Gaussian process Sh with 
covariance E[5^(t)5/ l (s)] = I(t, s). □ 

4. Properties of the limit process 

In this section, we study some properties of the limit process obtained in Theorem 2. First we deal with 
local self-similarity. 

Proposition 1. We assume that h is Holder continuous. Then the process Sh is locally self-similar, more 
precisely 



v 

lim 



i{ SUt + H; Shlt) } ={T( t ,«)},,„ >0 =T, (19) 

'I t J t,U>0 



and the tangent process T is the centered Gaussian field with covariance 

( R(h(t)Mt)) (U,\2h(t) I \„\2h(t) |„ „\2h{t)\ + 

E[T{t,u)T(s,v)) = \ 4h(t)--2h(t)U u \ y ' + M l« v\ y >) ift-s, 

Moreover, for every t, as e — > 0, the field {(Sh(t + eu) — Sh{t))/ £ h ^'~\u>o converges in distribution in the 
space of continuous functions endowed with the topology of the uniform convergence to T(t,-) which is a 
fractional Brownian motion, such that E[7~(i, l) 2 ] — R(h(t),h(t))/(2h(t) 2 — h(t)) with Hurst index h(t). 

Proof. We have 

-(S h (t + eu) - S h (t))(S h (s + ev) - S h (s)) 



E 



-h(t)+h{s) 
ft+eu i>s+ev 



I rt+eu ns-f-ev 

This) J, d6 J s toTVQ,*) 

1 r r 

J d0 j daTZ*(t + e6,s + sa) 

d0 dcrTZ*(t,s,e : a,e), 
Jo Jo 



£ h(t)+h(s) j t 
1 

£ h(t)+h 

fU pV 



where 



n*(t,s,e,a,e) = £h(t) l h[s) _ 2 n* (t+ s e,s + so-) . 

First, we assume that t = s. We have 

TZ# (t, t, 6, a, e) = e W+ee)+h(t+e,)-2h(t) 

x TZ(t + e6,t + so; h(t + e9),h(t + ea))\9 - a |M*+^)+Mt+^)-2, ( 2 0) 
Since TZ is continuous, when e — > we have for 8 =/= a 

TZ{t + s9,t + so, h{t + e0), h(t + ea))\9 - a \Ht+eO)+h{t+ecr)-2 

-^R(h(t),h(t))\6-a\ 2h ^- 2 . (21) 
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Moreover, h is Holder continuous, so there exist ah > and c/j > such that for every 9 £ [0, u] and a £ [0, v] 
\h(t + e9) -h(t)\<c h e ah u ah and |/i(t + ea) - h{t)\ < c h e ah v ah . (22) 

Hence, when e — > we have 

£ h(t+ee)+h(t+ea)-2h(t) = exp( (^ + £ #) + h( j. + £(t) _ 2 h(t)) log(e)) -> 1. (23) 

Combining (20), (21) and (23) we obtain that for every t and almost every 9 and a, when e^Owe have 

Tl*(t,t,6,a,e) -> R(h(t),h(t))\6 - al 2 ^" 2 . (24) 

Note that because (22) the convergence (23) is uniform in (0, er) £ [0,u] x [0,u]. Then using (20) and the 
fact that R is bounded, we prove that there exists a constant c such that for every e, t and almost every 
(6, a) £ [0,u] x [0,v] we have 

n*(t,t,9,cT,e)<cx{\9-a\ 2b - 2 + \9-<T\ 2a - 2 ). (25) 

Applying bounded convergence theorem and using (24) and (25), we get that the finite margins of {(Sh(t + 
eu) — Sh(t))/e h ^'} u >o converge to those of T{t, •) as e — > 0. It remains to prove the tightness. Let u < v such 
that | it — v\ < 1. We have 



E 



f S h (t + eu)-S h (t) _ S h (t + ev)-S h (t) Y 
\ £ h(t) £ h(t) J 

L-E[(S h (t + eu)-S h (t + ev)j 2 } = ^ ( - f + d9 f + daTZ*(9,a) 



\u — v\ 



2a 



£ 2h(t) 

<sup|i?| [ U d9 [ U d(r \e-a\^<-^^- 
Jv Jv 2a(2a - 1, 

Then by Kolmogorov criterium, the tightness holds. 
Now we assume that t ^ s. We have 

K*(t, s, 9, a, e) = £ ^- h (t)-h(s) n ^ + ^ g + £ ^ ^ + £ g^ h ( s + £(J ^ 

x\t-s + e(9 - (T) \Ht+ee)+h(s+e,)-2 ^ (26) 
Because ty^swe can check that for every 9 and a we have as e — > that 

K*(t,s,9,a,e)->0, (27) 

and for e sufficiently small, TZ^(t,s,6,a,e) is uniformly bounded for (9, a) £ [0,u] x [0,v]. Then by the 
bounded convergence theorem, we conclude the proof. □ 

It is classical to deduce pointwisc Holder continuity from local self-similarity [5]. 

Proposition 2. The process Sh admits a continuous modification. Moreover, if we assume that h is Holder 
continuous, then for every to £ M + the pointwise Holder exponent a~ (to) of Sh is almost surely equal to 
h(t Q ). 

Proof. We deduce as in [5] from Proposition 1 that a~ (to) < h(to). Now we prove that a-x (to) > h(to). 
We let < 77 < 1/2. For every s and t £ [to — ?/, to + n] such that s < t, we have 



E[(S h (t) - S h (s)f] < sup \R\ J d6 J da \9 - a\ 2inf ^ h ~ 2 

< sup \R\\t - s | 2inf M h < sup \R\\t - s | 2inf [*o-^*o+,i \ 
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By the fact that Sh is Gaussian and applying Kolmogorov theorem [7], we get that a— (to) > inf [ to _ I) to+); ] h 
for every < rj < 1/2. Then letting 77 — s- and using the continuity of h, we get a~ (to) > h(to), and hence 
a~ h (t )=h(t ). □ 

Please remark these properties are true even if (4) is not fulfilled. 



5. A representation of the limit process 

The aim of this subsection is to give a representation of the limit process Sh obtained from Theorem 2. This 
representation uses a universal Gaussian process that we introduce in the following section. 

5.1. An universal Gaussian field 

Here we consider a centered Gaussian field X = {X n (H)} nH . We define S N (t,H) for every N, t and H by 



1 L^vtJ 

S N (t,H) = - w J2 X ^ H )- 



(28) 



Theorem 3. Let R be an asymptotic covariance and {X n (H)} ri: H S Qr be a centered Gaussian field. Then as 
N goes to 00, the finite- dimensional distributions of {S N (t, H)} t ,H converge to those of a centered Gaussian 
field {W(t, H)}t,H with covariance given for every Hi, H2, t and s by 



E[W(t,H 1 )W(s,H 2 )} = 



R(H 2 , Hi) 



(Hi + H 2 )(Hi+H 2 -l) 



R(H U H 2 ) 



^ + H 2 )(Hi + H 2 - I) 

K{t,s;H u H 2 ) 
(H 1 + H 2 )(Hi+H 2 -l) 

Proof. We have for every Hi ,H 2 ,t, s, 

[Nt\ lNs\ 

1 



H1+H2 



t-s 



H 1 + H 2 



(29) 



E[S N (t,Hi)S N (s,H 2 )} = —— J2 E nXkmXjiH*)], 



3=1 h=l 

and we let 

K( Sl t) = [ d9 I dan(e,a;Hi,H 2 )\e - a\ Hl+H2 - 2 . 
Jo Jo 

Let 77 > 0. Using (12) and a< h<b, there exists an integer M = M„ such that for \j — k\ > M and every N, 



N 2 



K(j,k;H u H 2 ] 



j-k 



N 



H1+H2-2 



N 2 



J-k 



<E[X k (Hi)X 3 (H 2 )} 



<—TZ(j,k;H u H 2 ) 



j-k 



N 



H 1 +H 2 -2 



N 



V 

N 2 



H1+H2-2 



j-k 



N 



H 1 +H 2 -2 



(30) 
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Because of Lemma 2, 

K( S) t)-r,L(s,t)<limM V V E[X )b (F 1 )X J -(H 2 )]l b _ fe | >M 

3=1 fc=l 

[ivtj L^sj 

< limsup ^ ^ E[X fc (i7 1 )X J ( J ff 2 )]l b _ fc | >M < A-(«,t) + r/L(s,t), (31) 

AT — >oc -1,1 
J=l fc=l 

where 

L(s,t)= f d0 / da|e-a| ffl+H2 - 2 . (32) 
Jo Jo 

We can conclude as the end of the proof of Theorem 2, we get that 

lim E[S N (t,H 1 )S N (s,H 2 )} = [ dO [ da \0 - a\ Hl+H2 - 2 TZ(9, a: H u H 2 ). 
w ->°° Jo Jo 

Then by a direct computation of these last integrals, we get the convergence to a covariance given by (29). □ 

It is classical to give an alternative form of Theorem 1 which is based on a renormalization group (see, 
for instance, [15], pages 338-339). Now we propose this alternative approach for Theorem 3. 

We fix an asymptotic covariance R. Because of Theorem 3, there exists a process W = {W(t,H)}t,H, 
which is unique in distribution, such that for every {X n (H)} n _h £ Gr 



W 1 ^ 1 

I n=l ) t,H 



lim 

N 



We have the equality in distribution for every a > 

{W(at,H)} tiH = {a H W(t,H)} ttH . (33) 
We define the field Z = {Z n (H)} ny H for every n and H by 

Z n {H) = W(n, H) - W(n - 1, H). 

It can be verified that {Z n (H)} n) H & Gr- 

Now we define the renormalization semi-group Tn for every N. We let for every X G Gr 

T N X = {{T N XUH)} nH , (34) 

where for every n,H, 

(n+l)N 

(T N X) n (H) = ^ w ^ X n (H). (35) 

j=nN+l 

Because of (33), Z is a fixed point in Gr of T/v for every N. Moreover, for every X in Gr we have by 
Theorem 3 

lim {(T N X) n (H)} H = {Z n (H)} 

Finally, we have proved the following result. 

Theorem 4. The renormalization semi-group Tn admits the field Z as unique fixed point in the space Gr- 
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5.2. Representation theorem 



Here we consider an asymptotic covariancc R, a field X £ Qr and a continuous multifractional function 
/i:R+ — > [a, 6] C (1/2, 1). Recall that by Theorem 2, there exists a multifractional process SV such that 

K. n=l ) t 

and thanks to Theorem 3, there exists a field W such that 

W i ^ 1 

I n=l J t,iJ 

In this section, a representation theorem establishes the link between the process St and the field IV. 



Theorem 5. We assume that the function R is three times continuously differentiable and the function h 
is two times continuously differentiable. Then we have the following equality in distribution 



{Sh(t)} t>0 = {W(t,h(t))- / h'{9) 



I x y " Jo w dH 
where the right-hand side of (36) is always defined. 



de } , (36) 

H=h(6) 



Proof. We can deduce from the assumptions of the theorem and Kolmogorov's criterium that the sample 
paths of (t,H) i — ► W(t,H) are (almost surely and up to a modification) continuous with respect to the first 
variable and two times continuously differentiable with respect to the second variable. Hence, in particular, 
the right-hand side of (36) is always defined. 

Now we deal with (36). We consider the field {Z n (H)} n ,H defined in previous subsection for every n and 
H by Z n {H) = W(n, H) — Win — 1, H). On one hand, because of Theorem 2, we have 

t> ( lNti 7 (h N \) 

WEir ={&(*) W (37) 

l n=l ) t>0 

On the other hand, we shall prove that 



t>0 

which, combined with (37), proves (36). We have 



t>0 



INt] N \ ([Nt] \ 

£ = £ J^mn, C) - W{n - 1, 0) 

n=l ) t>0 I ™ =1 ) t>0 



t>0 



= {/f(*)--fa (*)} t >o. ( 39 ) 
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where for every t and N 



A' 



N 



m N 

At 



and 



LAtj 

h N (t) = E (w( ^.h?) -w( 



N 



N 



Almost surely, 



lim I»(t) = W(t,h(t)) 

N—>oo 



and thanks to Lemma 3 and regularity of W and h, 

limj?(t) = J* ti (0)^(0, h(0)) d9. 
Combined with (39), this proves (38) and concludes the proof. 



□ 



6. Examples 

In this section, we give some examples of multifractional processes Sh that we can obtain from Theorem 2 
as limits of sequences {S^}n for a multifractional function h: [0, oo) — * [a, 6] C (1/2, 1). In all this section, 
we assume that h is Holder continuous. 

6.1. Fractional white noise model 



Let us first consider the case of the multifractional Brownian motion, which at first motivated this article. 
In constrast to fractional Gaussian noise, which is the fixed point of the renormalization semi-group, it is 
not the case for the increments of multifractional Brownian motion. In this example, we investigate what is 
the limit for increments of multifractional Brownian motion. More precisely, let us consider for every t € K 



w H (t) 



1 



-itx 



l 



(40) 



C(H) 

where B is the Fourier transform of a real Gaussian measure B and the constant C{H) can be written as 

- | e -iz_1|2 



C(HY 



\2H+1 



■ dx ■ 



TX 



HT(2H) sin(7tff) ' 



Please note that for each H, Wh is a standard fractional Brownian motion. Moreover, if h : [0, oo) — > (1/2, 1) 
is a multifractional function, then t \— > Wh(t) (t) is a multifractional Brownian motion. We let 

X n {H) = W H {n + l)-W H {n). 

We compute the covariance between Xj(Hi) and Xk(H2) for every j, k, Hi and H2: 

C{{H l+ H 2 )/2)\ A uHi+H2 



K[Xj{H{)X k (H 2 



C{H X )C{H 2 



H 1 +H- 2 



+ 



1 - 



j - k 



-2 



(41) 
(42) 
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By a Taylor formula, we get as u — > 

|1 + u\ H ^ = 1 + (#! + # 2 )u + ^(-ffi + #2)(ffi + ^2 - 1)« 2 + 0(u 3 ), (43) 

where the is uniform in [Hi^H^) S [a,6] 2 because 1/2 < a < b < 1. Combining (41) and (43), we get that 
the asymptotic covariance R of {X n (H)} nj u can be written as 

Applying Theorem 2, we get that {S^)n converges to the process Sh with covariance 

E[S h (t)S h (s)}= [ d6 [ S da\9-a\ h ^ +h ^~ 2 (h(e) + h(a)) 
Jo Jo 

C((h{8) + h{a))/2f 



x (h(6) + h(a) - 1)- 



C(h(9))C(h(a)) 



Now we assume that h is continuously diffcrcntiable. Using Theorem 5, we can write Sh as the sum of a 
multifractional Brownian motion (defined by (2)) and a continuously diffcrcntiable process: 

U) = W Kt) (t)-f Q h'iB)^\ B=m 6B. (44) 

Moreover, we can obtain from (40) the harmonizable representation of the limit process Sh- 
f°° { (° itx - !) 

ft (e»»-i) ( log U 

Note that to establish (45) rigorously, we use the fact that the map / i— ► f(x)B(dx) is an isometry (in 

particular, \f(x)\ 2 dx = E[\ f(x)B(dx)\ 2 ] ), and the expression of the covariance of the derivative 
8W H /dH (seeTe.g., [1]). 

6.2. Fractional ARIMA model 

In the study of the previous example, we have obtained a symmetric asymptotic covariance R. In this section, 
we present a model for which this symmetric property is not satisfied. We also aim to generalize results in 
[14], where an invariance principle is established for a class of nonstationary processes with long memory. 
These processes are generalizations of FARIMA process. At the end of Section 1, the authors wonder if one 
can extend their results in a continuous time setup. Here we address this question in the Gaussian case. 

Let de (0,1/2). Wc consider the (Gaussian) fractional ARIMA model (FARIMA in short) defined for 
every n £ N by (see, for instance, [15]) 

d ^ T(d + l) 

^ = 1,-^^-1, (46) 

where {gi}i^z are i.i.d. standard normal random variables. Here we let X n {H) = <p" 1/2 . Wc first establish 
a lemma. 
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Lemma 1. As n — > oo, we have 



sup 

di,d 2 6[o-l/2,6-l/2] 

Proof. We have 



,1-*-* f r(rf 1 + n + i)r(rf 2 + i) 



(l + tif'-V^'du 



< Ji{d 1 ,d 2 ,n) +j 2 (d 1) d 2 ,ri) + j 3 (d 1 ,d 2i n), 



where 



7i(di,d2,n) 
72 (dii ^2, n) 
7 3 (di,d 2 ,n) 



,1— di— da 



2=1 



z 



d3 _Jr (d!+n + Q 
(n + l)! 



(Z + n) 



di-l 



di-1 / , \ d 2 -l 



(H-u)*- 1 «*'- 1 du 



Let us begin with 73. The function /d lj( j 2 : u l— * (1 + 1- " d2 1 i s decreasing on [0, 00) so we have 

/ fd 1 ,d 2 (u)du< -VVdLdJ - ) < / fd 1 ,d 2 (u)du. 
Ji/n n f~[ W Jo 



Hence, 



sup 73(di,d2,n) < 



(a - l/2)n a - 1 / 2 ra- 



ti. 



Now we deal with 72. Using Stirling formula, T(z) ^ z ^oo v27te z z z 1 / 2 , we have for every i] > that there 
exists M v depending only on 77, a and 6 such that for every Z > M,, 



r(di + Q 

Z! 



di-l_X-di 



<r)l ai - L e 



We get for every n > M,, 

00 

72(rfi,rf2,n)<n 1 - d ^ d ^ 



1=1 

00 



< 7771 



1 — di — d 2 



So, for every 77 > 0, limsup^^ sup dl d2 72(^1, d 2 ,ri) < -qS where S = sup dl d2 / °° fd u d a (u) du < 00. Hence, 
lim JWOO sup dl d2 72(^1, d 2 , n) =0. We use a similar argument for 71. □ 



Using Lemma 1, the function R is 

-1 



iss 
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and the limit process Sh has the covariance 

ft r s ,q_ \h{0)+h{a)-2 

E[S h (t)S h (s)] = I d0 I da 



o Jo T(h(9) - l/2)T(h(a) - 1/2) 

X / ( 1 + tt )h(»V<7)-3/a ti h(flA«r)-3/2 du 



JO 



|h(e)+h(er)-2 , 

r(2 - h{9) - h{o)) 



T{h(0) - l/2)T(h(a) - 1/2)' 
Appendix. Riemann sum convergence type lemma 

Lemma 2. Let t, s > ; G : [0, i] x [0, s] — > R &e a continuous function and M > 1 . We /lave 



1 



iim tp> y y g 





j — k 


\N N J 


N 



l|j-fe|>M 



3=1 fc=l 

where the integral in the right-hand side of (47) is always defined. 
Proof. We let for every j,k,N and (6, a) E X 



dO / do- G(0, a) \ 9- a 
o Jo 



h(8)+h(a)-2 



, (47) 



C 1- Jl] 


j — k 


\N N J 


N 



hZ+hf-2 



Notice that 



1 

"FT- 



im] INs] 

EE 1 

3=1 fc=l 





fc|>M 








j - k 


(1 


' NJ 


N 



h^+hf-2 



1 



JO 

hence, it suffices to prove 

\im j d6 j daG N { 



i-k\>Mi 



h(8)+h{<j)-2 



(48) 



We shall use a dominated convergence argument. Since G is continuous, then for almost every (6, a) G 

[0,t] X [0,8], 



Iim GAr(fl,o-)=G(ff,ff)|fl-a| 



h(6)+h(cr)-2 



(49) 



It remains to establish the boundedness of \Gn\ by an integrable function. Since G is bounded, then for 
every (0, a) £ (^, ^] x (V>lH 



|Gjv(0,(7)|<sup|G|l 





j - fc 


2a-2 




26-2^ 


( 


AT 


+ 


AT 
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and since \j — k\ > M > 1, 





j-k 




j - k 








N 




N 



hence, for almost every (0, a) € [0, f] x [0, s] 

\G N (e,a)\<svL V \G\(\e~a\ 2a - 2 + \9-a\ 2b - 2 ). 

The map (9, a)t-+\0- <r| 2a ~ 2 + \6- <r\ 2b ~ 2 is integrable on [0, t] x [0, s], then thanks to (49) and dominated 
convergence, we get (48), which proves (47). □ 

Lemma 3. Let a < b G M. and two functions f : R+ x [a, b] — > M. and h : M.+ — > [a, b] . We assume that h is 
two times continuously differentiable and f is two times continuously differentiable with respect to its second 
variable. Then for every t>0, we have 



[Nt] 



71=1 



N \N 



N 



N 



= / h'(9)d 2 f(9,h{9))d9. 



Proof. By the Taylor formula, for every n and N, there exists t Hi N € [(n — 1)/N,n/N] such that 
/ n — 1 / n\\ ( n — 1 /n — 1 



N ' ViV 



TV 



1 ,,/n-l\ „ Jn-l ,fn-l\\ 1 

= ^hr Whr' fc hr )) + w^ 



(50) 



where 



^ra,JV — h"(t n>N ) d 2 f 



Thanks to the regularity properties of /i and /, there exists C > such that for every n and N, |-E n ,iv| < C- 
Then using (50), 



LJVtj 



iV VJV 



N \ N 



[Nt] 



2 f\ -^rr-,h 



N 



N 



N 



< 



Ct 
TV' 



We conclude the proof by combining classical Riemann sums convergence with (51). 
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